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Abstract 

Knowledge of thermal diffusivity data for high-temperature materials (solids and liquids) 
is very important in analyzing a number of processes, among them solidification, crystal 
growth, and welding. However, reliable thermal diffusivity versus temperature data, par- 
ticularly those for high-temperature liquids, are still far from complete. The main mea- 
surement difficulties are due to the presence of convection and the requirement for a 
container. Fortunately, the availability of levitation techniques has made it possible to 
solve the containment problem. Based on the feasibility of the levitation technology, a 
new laser flash technique which is applicable to both levitated liquid and solid samples is 
being developed. At this point, the analysis for solid samples is near completion and 
highlights of the technique are presented here. 

The levitated solid sample which is assumed to be a sphere is subjected to a very short 
burst of high power radiant energy. The temperature of the irradiated surface area is ele- 
vated and a transient heat transfer process takes place within the sample. This contain- 
erless process is a two-dimensional unsteady heat conduction problem. Due to the non- 
linearity of the radiative plus convective boundary condition, an analytic solution cannot 
be obtained. Two options are available at this point. Firstly, the radiation boundary condi- 
tion can be linearized, which then accommodates a closed-form analytic solution. Com- 
parison of the analytic curves for the temperature rise at different points to the 
experimentally-measured values will then provide the thermal diffusivity values. Sec- 
ondly, one may set up an inverse conduction problem whereby experimentally obtained 
surface temperature history is used as the boundary conditions. The thermal diffusivity 
can then be elevated by minimizing the difference between the real heat flux boundary 
condition (radiation plus convection) and the measurements. 

Status of an experimental study directed at measuring the thermal diffusivity of high-tem- 
perature solid samples of pure Nickel and Inconel 718 superalloys are presented. Prelim- 
inary measurements showing surface temperature histories are discussed. 
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THE PHYSICAL PROPERTIES OF LIQUID METALS 


6.6. EXPERIMENTAL DATA FOR THE VISCOSITIES OF PURE LIQUID 

METALS 

It is very difficult to state definitely the accuracy of viscosity measurements for 
liquid metals. Errors of ± l to ± 20 per cent would seem to be a fair estimate 
with the exception of a few metals. There are not many well-established data 
for liquid metal viscosities. Viscosity data for pure liquid metals are listed in 
Table 6.3. Data for others currently available are given in Table 6.6. 

As shown in Fig. 6.26, experimental viscosity values for liquid aluminium 
exhibit very large discrepancies, and the values calculated from eqn (6.40) 



taken from Iida, T. and 
Guthrie, R. I. L. The 



Science Publications, 1988. 
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Fig. 6.26. Viscosity of liquid aluminium as a function of temperature, as determined by several 
workers; (1) Jones and Bartlett (1952-3); (2) Yao and Kondk (1952-3); (3) Yao (1956); (4) RothweU 
(1961-2); (5) Gebhardt and Detering (1959); (6) Gebhardt, Becker, and Dorner (1954); (7) Lihl. 
NachtigalL and Schwaiger (1968); predicted from eqn (6.40). 
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Pulse Diffusivity Method 
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FEATURES OF THE FLASH TECHNIQUE 


Original technique: 

Proposed by Parker et ai [l] 2 in 1961. 
i®* According to the Science Citation Index, more than 300 
scientific papers have referred to Parker et al . [l] 

The features of this technique were reviewed by Taylor and 
Maglic 121 in 1984: 

A literature survey in the early seventies showed that about 
75% of the thermal diffusivity data reported within that 
period were obtained using the flash method. 

Variety of materials including metals, alloys, ceramics, 
composites, etc. 

b® Thermal diffusivity values (10' 7 to 10' 3 m 2 /s) have been 
reported in the temperature range of 100 to 3300 K. 

*®* Small disk-shaped specimen (6—16 mm in diameter and 
less than 4 mm in thickness) 
i®* Half-time transients within 0.04 to 0.25 seconds 

Variations of the flash method up to now account for: 

Radiation Losses Multi-Layer Specimen 

Finite Pulse Time Large Sample Effects 


1. W. J. Parker, R. J. Jenkins, C. P. Butler, and G. L. Abbott, "Flash Method of 
Determining Thermal Diffusivity, Heat Capacity, and Thermal Conductivity," J. App l. 
Phvs. . 1961, Vol. 32, pp. 1679-1684. 

2. R. E. Taylor and K. D. Maglic, "Pulse Method for Thermal Diffusivity Measurement," 
Compendium of Thermophvsical Property Measurement Methods . Vol. 1, edited by K. 
D. Maglic, Plenum Publishing Corporation, 1984, pp. 305-336. 
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THEORETICAL ANALYSIS 
Governing Equation 


Assuming that the spherical sample is homogeneous and the thermophysical properties 
are independent of temperature for the given temperature range of interest, the extended 
flash technique can be modeled as a two-dimensional (azimuthal symmetry) unsteady heat 
conduction within the sphere. 


j__8_/ r2 ae\ 

r 2 dr\ drj 


+ 


j__a_ 

r 2 d\i 


(1-H 2 )-^ 


80 

d\i 


= is®. 

a dt 


where: 

r radial coordinate within the sphere, m 

t time, s 

T temperature, K 

temperature of the surroundings, K 
oc thermal diffusivity, m 2 /s 

0 temperature difference, defined as (T - TJ, °C 
H transformed 0 coordinate, defined as cos 6 
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THEORETICAL ANALYSIS 
Assumptions Related to the Initial Conditions 


The following assumptions are made: 

1 . The intensity distribution of the radiation energy flux of the incident beam is known. 
In cases where a laser is utilized, a Gaussian distribution is assumed: 

Q(r b ) = Q m exp(-2r b /w b ) 

where: 

Q m peak energy flux at the center of the beam, J/m 2 
r b radial distance from the center of the beam, m 
w b diameter of the laser beam, m 

2. The energy radiation is absorbed instantaneously and uniformly within a thin layer 
with the thickness g at the irradiated front surface (r = R, 0 <, 6 ^ 0 O , 0 £ <f> < 
2t). This layer of material is very thin compared to the radius of the sample. 

3. The duration of the radiation burst is negligible compared to the rise time to half- 
maximum temperature characteristics being sought. 


THEORETICAL ANALYSIS 
Initial Conditions 


The initial condition can be written as: 

8(r*,p,0) = — ^ (( 1 -J?) z r* z 1, p 0 s \i z 1) 

P C p g R ror 

6 (r*, ii,0) =0 rest of the sphere. 


where the assumption T ( = T^ is employed. The angle d Q represents the relative size 
of the sample and beam diameters. 


NOTE: More realistic initial conditions are currently being studied. 
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THEORETICAL ANALYSIS 
Boundary Conditions 

The appropriate boundary condition accounting for both radiative 
and convective heat losses is: 

-k T„ = € eff a F (T 4 - o + h (T - T J. 

Upon linearization, it can be expressed as: 

{R, p, t) = (4 e eff aFTsur + h) 0(i?,p,t) 

where: 

k thermal conductivity, W/(mK) 

e eff effective hemispherical emissivity of the sample 

F radiation view factor 

a Stefan-Boltzmann constant 

h convection heat transfer coefficient. 


Alternately, by defining r* = r/R, we can write: 

( 1 , 11 , t) = Bi 0(1,11, t) 

dr 

The Biot number (Bi) which is defined as: 

_ . _ (4£ e ff aFT sur + h) R 

Bl = k 

is a measure of both convective and radiative surface heat 
losses. In the case of no heat loss, the Biot number equals zero. 
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THERMAL DIFFUSIVITY DETERMINATION 


A simple one-step method to determine the thermal diffusivity 
which is independent of the heat loss parameter (Bi) is proposed. 

Record the temperature rise histories, at least at two points 
on the surface of the sphere simultaneously. 

Determine the rise time to the half-maximum temperature 
from the experimental temperature rise curves. 

^ For a given energy distribution and incidence angle, the 
relationships between and Bi at the two points are 
expressed as: 

[ co i/2 J i = a ± e ^ lv/ST at point P 1 


2 



e ~b 2 VBi 


at point P 2 


Because the Biot number should be the same in both equations, 
it can be eliminated. Finally: 


a 




[* 1 / 2 ]^ 


where the parameters A and B are as follows: 


A = 



B = 


b 2~bi 
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The Values of A and B 
(?! = 90° and P 2 = 180°) 


Gaussian Energy Flux Distribution 
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Uniform Energy Flux Distribution 
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Figure 2. Photograph of experimental apparatus. Left to 

right: data acquisition computer, line-scan 

radiometer, hot specimen in vertical glass tube 
with laser window and gas supply, optical 
pyrometer, induction heater power supply, manometer 
and gas handling equipment. 
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12 mm Inconel Specimen - File in12.10 



Time, S 


6 mm Nickel Specimen - File ni06.26 




Time, S 


Figure 5. Temperature-tine response for pulsed heating 

experiments on a 1.2 cm Inconel and a 0.6 cm 
diameter nickel sphere. The plots show the 
normalised temperature determined from the 
calibration pixels as a function. 
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CONCLUSIONS 


«r In order to determine the thermal diffusivity of high-temperature materials, an 
extended single-step containerless flash technique which is applicable to levitated 
spherical specimen is proposed. 

» The thermal diffusivity is determined by knowing the sample diameter and recording 
the temperature rise at least at two different points on surface of the sphere 
simultaneously. 

w The main advantage of this method is that the thermal diffusivity can be determined 
without the necessity of knowing heat losses at the surface. In addition, the limiting 
stringent specimen size requirements associated with similar methods is no longer 
applicable. 


NEXT PHASE 

Extend the theoretical linearized analysis to other liquid droplet shapes (e.g. ellipsoid 
of revolution) 

«r Solve the "Inverse Problem" to determine thermal diffusivity 
**• Explore Marangoni (surface-tension-driven) convection effects 
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